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We formulate the dynamical theory of nematic polymers, starting from a microscopic 

O I Poisson bracket approach. We find that the Poisson bracket between the nematic director 

and momentum depends on the (Maier-Saupe) order parameter of the nematic phase. We 

►^ I use this to derive reactive couplings of the nematic director to the strain rates. Additionally, 

2^ . we find that local dynamics breaks down as the polymers begin to overlap. We offer a 

^ ; physical picture for both results. 

O 

a^ 
a^ 

-(— > 

a 

o 
o 

> 

X 



21 June 1999; revised 29 June 1999 
t email: kamien@physics.upenn.edu 



1. Introduction and Summary 

The dynamics of line-lilce objects, from biomolecules to fiux-lines in superconductors, 
has become increasingly important to understand. The processing and manufacture of 
petroleum products, food products and petroleum-food products requires the control of 
polymer flow. The remarkable properties of superconductors can be exploited only if flux- 
flow dynamics is controlled. While the equilibrium statistical mechanics of line-like objects 
such as directed polymers, flux lines and dipole chains in ferro- and electrorheological fluids 
has been formulated from the microscopic physics [jl|,|^, a similar theory has not been de- 
veloped for the dynamics. While polymer entanglement dramatically affects the dynamics 
by introducing cosmologically long time scales, the statics is not affected. Indeed, it is the 
dynamics of entanglement that aids in pinning of high-temperature superconductors, in the 
strength of glue and in other remarkable rheological properties. Previously, dynamics have 
been formulated on a hydrodynamic basis, which ignores the connectivity of the lines ||^ . 
Other approaches include effective theories which use single-polymer response functions 
and neglect polymer-polymer interactions [Q. Here, we study the dynamics of polymers 
which, due to steric or other interactions, are already aligned into a nematic state. By 
considering a system which is orient ationally ordered in equilibrium, the effect of flow on 
alignment can be separated from the effect of alignment on flow. 

In this work we formulate the fluctuating hydrodynamics of polymer nematics based 
on a microscopic, Poisson bracket approach [^. In this way we will account for polymer 
degrees of freedom as well as long-lived hydrodynamic variables. Our principle result is a 
fundamentally new derivation of the elusive reactive coupling A between the director 
fleld n and the velocity fleld v: 

-KT = i^f^'^ - n^n^] I ——n^d-yV^ —n^d^v^ [ . (1.1) 

Unlike most parameters in hydrodynamics A is neither set by Kubo formulae nor the 
fluctuation-dissipation theorem. We will show that to the order we work, our result 
agrees with the analysis of Archer and Larson [[7| which is based upon a phenomeno- 
logical. Landau-theory dynamics ^ . Thus our result both verifles the validity of both a 
microscopic and a phenomenological approach. 

As is typical of the former approach, we develop our theory in steps. First we identify 
the hydrodynamic variables. In a polymer nematic system we have four flelds of interest: 
the areal polymer density /o(x), the momentum density g(x), the fluctuation of the nematic 
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director 5n(x) = n(x) — no, and the density of polymer "heads" and "tails", Pht(x). We 
then write these fields in terms of the positions and momenta of the individual monomers. 
We consider, as an example, the areal density field p(x). Since the polymers are directed 
we may write the location of the a^^ polymer in terms of an affine parameter Tq, which 
marks the monomers: 

Ra(ra) = K{Z) + Tc,h''^,rl{z) + Tc,hl,Tc, + Sc,] (1.2) 

where fa{z) is the (two-dimensional) displacement of the monomer at height z = r^t + s^ 
away from the ground state, straight-line configuration with the two-dimensional tangent 
/la, and the polymers start at a height z = Sa- The expression for the areal density is: 

p(x) = ^5^ l^fa(z) -F(2:-Sa)/ia -x_L 6 [^ - s^] 6 [c^ - 2;] , (1.3) 

where Cq, is the height of the polymer end. This expression requires further explanation 
and, at the same time, shows the inherent complexity of this problem. The areal density 
is a sum of delta-functions in the plane at the positions of the polymers. Moreover, since 
these polymers have heads and tails the product of Heaviside step functions O(-) counts 
the polymer only if we consider heights between the start and end of the polymer at z = s^ 
and z = Cct respectively. We will investigate similar expressions for the remaining fields in 
the next section. 

We note here that although expressions such as ( |1.2| ) are awkward, similar expressions 
could not easily be formed in isotropic polymer systems. The essential feature of polymer 
nematics is that chain contour length corresponds to the common z-axis along which all 
the polymers are more or less aligned. This allows us to formulate our theory in space, 
so as to incorporate hydrodynamics, while at the same time including the microscopic 
polymer degrees of freedom. These degrees of freedom are accounted for via the Poisson 
brackets of the r'a{z) with the conjugate momentum Paiz). From these, we can generate 
the Poisson brackets of the hydrodynamic fields: if we can get a closed set of brackets 
then we can construct a consistent hydrodynamic theory. However, we see that the areal 
density ( |1.3D not only depends on fa{z) but also on the location of the polymer heads 
and tails as well as the average polymer tangent /Iq. The head and tail locations are 
accounted for via the field Pht(x) described above. What do we do with the average poly- 
mer tangent /iq,? Though the nematic phase has two Nambu-Goldstone modes associated 
with spontaneously-broken rotational invariance, there is, in addition, a massive mode, the 



Maier-Saupe order parameter 5" which measures the degree of nematic ordering. The ha 
are associated with this mode. The order parameter can change locaUy via changes in the 
local polymer alignment. Presumably, if the polymers are entangled, the time scale for 
these rearrangements corresponds to some sort of reptation time required for the release 
of topological constraints. On the other hand, if the polymers are short (unentangled) 
then there is a separation of time scales which allows us to approximately decouple the 
local Maier-Saupe mode. In both these limits, we can decouple these tilt modes from the 
monomer modes. This will allow us to pre-average the tilt modes to arrive at effective 
Poisson brackets for the remaining fields. 

The next step in our formulation of dynamics requires a coarse-grained free energy 
in terms of our macroscopic variables. In section III we will propose a free energy con- 
sistent with the symmetries of the system which also includes the equilibrium physics of 
the polymer ends. From this we will derive dynamical equations for the hydro dynamic 
variables. 

2. Definition of Hydrodynamic Variables 

We start by formulating the theory in terms of microscopic, directed trajectories of 
the individual polymers, {fa{z)}, where a labels the polymer and r{z) is the displacement 
of the polymer from its equilibrium position. Recall that the nematic state is characterized 
by the Maier-Saupe order parameter, 

s^^A^^ill^ (2.1) 

where (■) denotes a thermal, ensemble average, and Q is the angle between the nematic 
director and the ordering axis (taken throughout to be z). When S ^ \, the directors 
do not all line up along z. Since the nematic order parameter is set by the details of the 
isotropic-nematic transition, it is a massive, non-hydrodynamic mode. In order to take into 
account the ground state value of S we will assume that the equilibrium polymer trajectory 
is Rq(2) = [(2 — Sq)/i^, (z — Sq)/i^,2] and that (rQ,(z) ) = 0. Thus, the microscopic director 
is dr/dz -}- h. Finally, the /U component of the momentum of polymer a at ^ is p^ [z) . 

Unlike previous formulations , we do not consider the traceless symmetric tensor Qij 
typically used for nematics. Rather, deep in the nematic state, we are only concentrating 
on the two Goldstone modes of the broken rotational symmetry. As we shall see this leads 
to a closed set of Poisson brackets to lowest order in derivatives and the director fluctuation 
6n. We will incorporate the "up-down" symmetry of nematics via the symmetry of the 
free energy, just as is done in the usual equilibrium description of nematics in terms of the 
unit- vector director field fi. 



2.1. Canonical Variables and Separation of Time Scales 

While our description of the chain locations is adequate in terms of r^^z) and /Iq, we 
should take care in finding correct canonical variables from which to construct macroscopic 
Poisson brackets. We imagine starting with a Hamiltonian which is a function the actual 
polymer location RQ,(rQ,), as in (|L^), and the conjugate momentum Pa(rQ,). The equations 
of motion in the xy-plane are 

If hoi is constant then these equations of motion are precisely 






(2.3) 



with 

KW,r^(^')} = {PaW,PM^')} = 0' (2-4) 

and 

K(^),pJ(/)} = 5„^5'^5(^-/), (2.5) 

where i and j run over the xy-plane. We have made the identification Pa(^) = Vaija) 
with z = Tq + Sq. We must also calculate the Poisson bracket of p^(z) with the spatial 
coordinate. Motion along the 2;-axis should be accounted for in motion of the polymer 
ends. Indeed, motion of the monomers up or down depends only on the motion of the 
starting point s^- Thus 

K(z),p^(/)} = (2.6) 



but 



{zc,p^p{z')] = 5cp5{zoi- z), (2.7) 



where we use the symbol z^ to remind the reader that the commutator does not vanish 
only for momenta and positions on the same polymer. In other words, z^ takes the place 
of the z-component of ri{z) and Za = Sa + Ta. We remind the reader that the fact that 
the affine parameter r can be mapped one-to-one with the coordinate height z is precisely 
the simplification which we are exploiting. Of course, when we consider hydrodynamic 
variables, z^ will just become the height coordinate z. 



We must now justify our neglect of the dynamics of ha- As we discussed in the 
introduction, we can do this in two Hniits. One hmit is when the chains are long and 
highly-entangled. In this case we expect that rearrangements of the nematic texture require 
topological constraint release or, in other words, reptation. We are concentrating here on 
dynamics on time scales shorter than the reptation time. We could extend this work by 
including some sort of reptation-based dynamics for the tilt fields ha, thus incorporating 
both hydrodynamics and topological constraints in the same theory. 

There is another limit of interest. If the polymers are short, we expect that the 
modes associated with the Maier-Saupe order parameter will be the slowest to relax to an 
equilibrium distribution. In this case, we expect that on the intermediate time scales we are 
focusing on, the distribution of ha will also be static, but consistent with the background 
nematic order. In this case the tilt modes can be averaged over, thus explicitly removing 
the dependence of the fields on them. We contrast this average to a thermodynamic 
average: in this case, the average is over an ensemble of different frozen-in textures while 
in thermodynamics we take the time average of a single molecule. In both these cases we 
get the same mean-square average of the tilt field, though the correlations will differ. 

Thus, in both the very long polymer and short polymer case we can pre-average the 
tilt degrees of freedom. 

2.2. Macroscopic Fields 

We must first define the coarse-grained, macroscopic fields, namely the mass density, 
the director and the momentum density. The mass density is given by: 

p(x) = ^3"^ [Taiz) + iz- Sa)ha - X_L j Q [z - Sa] 6 [Ca " z] (2.8a) 

As described in the introduction, the terms in the sum come from the following considera- 
tions: 1) The Dirac delta function "counts" the areal mass (note that it is a two-dimensional 
delta function, defined in the xy plane). 2) The pair of Heaviside step functions (0[-]) con- 
trol the polymer lengths. The a^^ polymer begins at ^ = Sa and ends aX z — Ca- 
The remaining macroscopic fields are defined via 

[p5n^](x) = J2 f ^ + ha) S' (fa{z) + {z- sa)ha - xA ^l^^-^^lSfc^^ (2.86) 
a ^ ^^ ^ ' ^Jl + hl 

(7^(X) = 5^K(^)5' (faiz) + ha{z - Sa) " X±) Q [z - Sa] 6 [Ca ' z] (2.8c) 



/ ^ \-l/2 

Note that an additional factor of ( 1 + /i^ 1 arises in the definition of 5n(x) since the 

normahzed unit (average) tangent vector is 

T, = I^^^^^M. (2.9) 

Note that the naive definitions of the macroscopic fields, e.g. 

p(x) = 5^52[f, - x^]e [z - s„] e [e« - z] , (2.10) 

a 

applies only when h^ = 0, or, in other words, when the nematic ground state has Maier- 
Saupe order parameter, S = 1. 

3. Pre-Averaging the Nematic Texture 

As we discussed in the introduction, the Poisson brackets of the hydrodynamic fields 
cannot close since they all depend on the average polymer tilt /iq,. As discussed in the 
preceding, we would like to pre-average over these degrees of freedom to arrive at Poisson 
brackets for the hydrodynamic fields. This procedure can be justified in two limiting cases. 

Consider the dynamics of a single polymer of length L. This can be decomposed into 
normal modes aia Rouse or Zimm [§ . In either case the n^^ mode has wavenumber q^ = 
2nn/L with a relaxation frequency u^ oc n'^/L'^. For very short polymers the separation in 
time scales between the lowest mode n = 1 and the other modes is large. In this limit we 
can decouple the longest mode - the tilt mode - characterized by h. We could thus average 
the tilts over the collection of polymers. The mean-square of the tilt fields is simply related 
to the Maier-Saupe order parameter S. We make the distinction between the time average 
of a single polymer's tilt field which gives equilibrium statistics and the quenched average 
over the frozen polymer tilts. 

As we imagine lengthening the polymers, the splitting in timescales becomes less pro- 
nounced and their spectrum becomes continuous. However, when the polymer is sufficiently 
long, a new dynamics dominates the longest wavenumber modes: reptation dynamics. In 
this regime, the tilt modes can only relax via activated reptation |]l^ which is an especially 



apt description of the likely rearrangements in polymer nematics. Thus, when the polymer 
is long enough to be entangled we also have a large separation of time scales. As a result, 
in this limit, we may also treat the polymer tilt field as essentially frozen. 
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We have now two distinct averages: the first is an average over a massive degree 
of freedom that sets the value of S. In the present formulation, this average is over 
hot- The second average is the thermal average over fluctuations about the ground state, 
i.e. averages over ra{z). We will treat these two averages separately, which amounts to 
treating the average over ha as a quenched average. To connect with the Maier-Saupe 
order parameter, we have 

since cos 9 in (|2.1|) is just the z component of the average tangent vector and where we have 
denoted the quenched average of X by X. Expanding for small \ha\, we have |(1 — S") ~ /i^. 
We will then take averages over h^ weighted by 

PiK) = ^^-"^''/'^ (3-2) 

where A = -1(1 — 5") so that (|3.1| ) is satisfied. 

In order to find the average macroscopic field variables, we must regularize the Dirac 
delta functions appearing in (|2.8|) . To do this we represent them as Gaussians with the 



widths taken to be the excluded area of each polymer in a fixed 2;-plane. We define 

52(r;a) = -^e-|^"*lV2a. (3.3) 

zna 



Then, when we average over ha weighting by ( p.2|) we find 



5^[fa{z) +{Z- Sa)ha - X_l; o] = 5^[fa{z) - X^; tt + A(z - Sa)\ (3.4) 

We can thus calculate the tilt-averaged fields. Since we have only included terms for 
the tangent vector in ( p.86|) to leading order in ha, we calculate average field values to 
order A: 

p(x) = ^ 5^ (^a(^) - x_l; a + A(2 - Saf) Sa(2) (3.5a) 

a 

[p5n^](x) = 5^ (1 - A) -^S' {fa{z) -xr,a + A{z- Sa)') S«(^) 

a 

- Y,i^ - Sa)Aa,52 {ra{z) ~xr,a + A{z - Saf) ^c.{z) (3.56) 

a 

ri^) = $^<(^)5' {rM -xr,a + Aiz - Sa)^) ^M (3.5c) 

a 

where Sq.(z) = O [z — Sq.] O [cq. — z]. 



3.1. The Breakdown of Local Dynamics 

Typically, Poisson brackets for dynamical variables are of the form 

{p(x), g{x')} = -V, [p{x)5\x - x')] . (3.6) 

The presence of the delta-function in (|3.6|) along with a local free energy density leads to 



local dynamics [T^], i.e. every term in the dynamical equations is evaluated at the same 
point in space and time. Such a result must, of course, be taken with a grain of salt. 
The delta function is an idealization to a system composed of point-like constituents. In 
reality, these delta functions should be replaced with a smeared-out distribution. Usually, 
this presents no problem - calculations may be done with the smeared-out distribution 
and the limit may be taken at the end. However, this is not the case in this situation. 

From (|3.4|) , it is easy to see that we may only replace the disorder smeared delta- 
functions with the original distributions (of width a) when a ^ ^{z — Sa)^. In other 
words, it is only when the original width is wider than the disorder-averaged width that 
the distributions are delta-function like. Since the presence of the term G [z — Sa\ O [ca — z] 
in the expressions for the hydro dynamic variables limits the maximum value oi z, we have 
^{z — SaY < ^(ca — SaY = A£^, where I is the typical polymer length. Thus we see that 
the original distributions are recovered everywhere whenever £v /i^ = ^vA <^ a. This is 
easy to interpret; when the average wandering of the polymer away from its starting point, 
£vA is smaller than the interpolymer spacing, we may continue to treat the polymers and 
their interactions as local. Once the polymers have tipped out of their cages the dynamics 
will become, necessarily, non-local. A force at x leads to a reaction at x' if a single polymer 
can go from x to x'. Presumably, the typical polymer length i should be replaced by ip, 
the polymer persistence length, when i > ip and the locality condition becomes A£p -C a^. 

In this delocalized regime the formalism will break down. In particular, the Pois- 



son brackets of the coarse-grained fields, calculated in terms of ( p.6|) and (2/7) will be 
straightforward but complicated, and will necessarily lead to non-local dynamics, even if 
the free-energy density arises from local interactions. In the Conclusion we will discuss 
possible alternatives to a non-local formalism based on different sets of variables. 



4. Fluctuating Hydrodynamics 

We can, however, consider the case of A ^ 1. This could arise in the description 
of flux lines in superconductors or polymer nematics in applied fields where the ground 
state polymer configurations are almost always parallel to z. In any of these limits, it is 
appropriate to replace the spread out delta functions 5^ (•; a + A(2 — s)^) with point-like 
delta functions. In this case, we can calculate the Poisson brackets of the average fields in 
terms of the canonical, microscopic brackets. We find the following, non-zero brackets, to 
lowest order in derivatives (of both Vfy^i^z) and delta functions) and 5n: 



{p(x), ^^(x')} = - [5^ - b^Jl\ dl {p{^)8\x-x') 

+ (1 + A) 5^ [pHT(x)53(f - x') - d, {p5n'{^)5^{x - x'))] (4.1a) 
{5n^(x),^^(x')} = -[K- d^Jl] 5\x - x')d^Jn\^) 

+ [(1 - A) 5'^dl - /^5^'dl] 5^{x - x') (4.16) 

{^^(x),r(x')} = -a: [g^{^)5%x-x')\+d^^, [g'' {^)d%x - x')] (4.1c) 

We have been forced to introduce an additional hydrodynamic field, Pht(x), the density 
of polymer heads and tails, defined as 

Pht(x) = ^ (1 - A) {52(fa(sa) - xi_)5{z - Sa) - 5^ra(ea) - XA_)5{ec, - z)] (4.2) 

a 

Note that the second term in ([4.1a|) is exactly equal to dzp{^)- In other words, we have 



(1 - A) a,p(x) + a, [p5n,(x)] = Pht(x) (4.3) 

This constraint is the familiar constraint of line liquids [P!^,p!^ which conserves polymer 



number but for polymer ends. The prefactor of (1 — A) comes from the average over the z 
direction. In other words, the rotationally invariant conservation law is no-'Vp+V-n = Pht- 
When averaging over the intrinsic randomness of no, we get a factor of (1 — A) in front 
of dzP- It is not present in other terms because they are higher order in derivatives and 
powers of 6n. 
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4.1. Equations of Motion 

With these Poisson brackets we can systematically derive the reactive terms of hy- 
drodynamics. The viscosities will appear in accord with the symmetries of the system. 



Being uniaxial there will be 5 viscosities [jlj], characterized by the viscous stress tensor 
Va/3 = ?7a/37/3^7/3 where 

where 5^^, = d^i, — n^rijj. The form of V is fixed by symmetries: by nematic inversion 
it must be even in n and by the Kubo formulae § it must be symmetric under fj, ^->- u, 
7 ^> /3 and (/xz^) ^^ (7/3), and A^p = ^ {d^vp + dpVj) is the strain rate tensor constructed 
out of the fluid velocity v. Expanding n = z + 6fi we can flnd the viscosity tensor in the 
nematic phase. 

The free energy from which we will derive the reactive terms is: 

F' = jd\ ||_ + Feq,ii| (4.5) 

where the equilibrium free energy is the sum of terms: 

-^equil = FsH + Fpol (4.6) 



and 



3 J Ki 2 , K2 , „^s2 . Ks 2 



Fsn = I d'x {^ {V-6nr + ^ (V X SnY + -f {dMY ) (4.7) 



is the usual Frank free energy for a nematic, while 

i^poi = I d\ ||5p2 + ^ [d,5p + poV-5nf\ (4.8) 

is the free energy of the polymer p|. In (|4.8| ) B is the two-dimensional bulk modulus, 
Pq = p — 5 p \s the average, areal polymer density and G is the fugacity for hairpins and 
free ends. We have replaced Pht with the expression in (|4.3|) , to lowest order in the fleld 
fluctuations 5p and 5n (which is why we drop the factor of (1 — A) in front of dzSp). The 
parameter G is given by G = k^Tl/2pQ where I is the typical polymer length [§. Terms 
of the form 5pV -dn are not allowed due to the n —>■ —n symmetry of the nematic phase. 
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We are now able to derive the equations of motion for g, 5n and 5p (where we have 
dropped the bar over the variables denoting the h average). They are (with v = g/p): 

dtp = -digi + (1 + A) Pht^^z - di {Srug^) (4.9a) 

SFlSn] 
dtSn, = -dj ivj6n,) + (1 - A) d,v, - Ad,v, - T— ^ + 6, (4.96) 

d{dni) 

dtg^ = -d, (^^ + [6i: - 6^6^] pd, (^) + d^V^^ + i^ (4.9c) 

where, in order to approach thermodynamic equilibrium, the noise terms ^i{x,t) and 
9i{x,t) have the following correlations ||15|| : 



( U^. t)^p{x', t') ) = 2k^Tr]^,^pd,d^5^{x - x')S{t - t') (4.10a) 

( e^{x, t)e^{xt') ) = 2k^TT6-l;J^{x - x)S{t - t') (4.106) 

{^^ix,t)e,{x',t')) = (4.10c) 



There are a number of features of these equations worth mention. First, we see that 
the typical polymer lengths come into the equations through the coupling G, present in 
the equations for dtSui and dtg^. When G is large, corresponding to long polymers, we 
may take the coefficient, Pht = 0. In the short molecule limit, presumably G ~ 0. There 
is a crossover wavevector kc = ^j2Bf^/{k^/Tt} below which the system behaves as a pure 
nematic and above which the length of the polymers becomes important. There is a 
corresponding timescale cUc oc /c^ ~ l/£ which delineates a similar short-to-long crossover. 

In the short-polymer limit, care must be taken due to the nature of the approximation 
under which the Poisson brackets were derived. In particular, since we treat the polymer 
tangent field as a vector, we are, in principle, distinguishing n and —n. However, for small 
fluctuations around equilibrium, long polymer nematics should not sample the entire space 
of fluctuations: it is very unlikely that an entire polymer will rotate by tt around the x- 
or y-axis. This simplification is what enabled us to get a closed set of Poisson brackets 
order by order in a derivative expansion. The effect of hairpins could be put in explicitly 
through the introduction of a hairpin- density field . 
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4-2. Reactive Couplings in the Nematodynamic Limit 



Finally, we note that (|4.96| ) gives an expression for the elusive reactive parameter A 



1^ in nematodynamics. By rotational invariance, n — » — n and the constraint n^ = 1, the 
reactive part of the equation of motion for n has the form: 



dt ~ ^"' 
We can thus identify 



1 + A _ill^ r. 



(4.11) 



A = i-2A = ^^^. (4.12) 

This is in general agreement with the results of Forster [^ who calculated the value of A 
within a Poisson bracket formalism for Q^i, the symmetric, traceless order parameter for 
short molecule nematic liquid crystals [0. He found that 

A ^ ii|^, (4.13) 

where a = {Ii + 2It)/{Ii — It) is a parameter depending on the moments of inertia /; and 
It of the nematogens parallel and perpendicular to the nematic axis, respectively. We thus 
recover Forster 's result when ct = 1 or, in other words, when the aspect ratio Ii/It of the 
nematogens becomes infinite. Presumably this is a consequence of taking delta-function 
densities in the transverse plane. 

It is interesting to compare this result to the work of Archer and Larson 0. In 
the same limit of infinite aspect ratio, they found an expression for A in terms of the 
expectations of the second- and fourth-rank order parameters, {P2) = S and ( P4 ) where 
P2 and P4 are the second and fourth Legendre polynomials evaluated at a; = cos 9: 

^^15(P.)+48(A)+42_ 

105(P2) ^ ^ 



To be consistent, we should expand (|4.14| ) in powers of A = (1 — S')/3 to compare with 



( [4.12| ) and compare linear terms. First, to linear order in A, 



;cos^^) = l-2{hl) ^ (l-ihDy ^ {cos^O)^ (4.15) 



and so ( P4 ) ~ (35^^ _ iQg _ y^/^g ^i^j^i^ g^^^gg 

2,„ .. 2 



X = 1 + Z^S-1) + ^{S-1)' + .... (4.16) 
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Thus to leading order in (S" — 1) = 3 A ( 4.14 ) agrees with (|4.12 )! Thus within the 5" ~ 1 



hmit our result should be consistent with the data as in 0]. We note, moreover, that 
since the isotropic-to-nematic phase transition is first order, S does not grow continuously 
from 0. Indeed, in Maier-Saupe theory S ~ 0.44 at this transition |lj], and thus A ^ 0.2. 



Therefore the deviation between our result and the more exact result (|4.14|) should be 
small sufficiently well aligned samples. However, although (|4.12| ) may be quantitatively 
reasonable, it misses an essential qualitative feature: it is always less than 1 and thus 
predicts that nematics will always tumble. 

Though the linear result in A fails to predict the crossover from flow aligning to 
tumbling behavior, the virtue of our derivation of A is that we get a direct interpretation 
of its origin. In a highly aligned sample with S" = 1 only gradients of v± along the z- 
direction can lead to rotations of the molecule (see Figure 1). On the other hand, when 
S < 1, gradients in the x-direction of Vz can also rotate the molecules and so the polymer 
nematics would always be in a tumbling mode. A higher order analysis in powers of A 
would be required to see if a A > 1 could come from our direct approach. 

5. Conclusions 

In summary we have derived the Poisson brackets of the relevant degrees of freedom 
for a polymer nematic. We have shown that these canonical brackets become highly non- 
local when the polymers start to overlap. In the limit where the polymers do not overlap, 
we have presented a microscopic derivation of the reactive coupling A, a coupling which is 
not a long- wavelength limit of a correlation function. 

Our analysis requires a "locality condition" in order to be trustworthy. While this may 
appear restrictive, one might imagine that, at some level of coarse graining the polymers, 
the no overlap condition might be met. If we were to clump polymer regions correlated 
in the xy-plane together and coarse grain to the scale of this "entanglement correlation 
length" ^e, our analysis may be applicable. Our analysis is especially applicable to flux 
lines in superconductors. There, Pht = and there is no problem, even in principle, to 
considering the flux line tangents as vectors. The flux lines, in addition, have A = 0, 
simplifying the theory further. It is perhaps in this context that one could hope to make 
the most theoretical progress. 

This derivation has shown the inherent, unavoidable no n- locality in polymer nematic 
dynamics. An interesting possibility is to focus on a different set of conserved or almost 
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conserved variables. For instance, it may be possible to reformulate this dynamics in 
terms of an entanglement density by calculating the Hopf density [jl^, a scalar which 



measures the local curvature of the director configuration - for long polymers curvature 
is a measure of local entanglement. This density could be used as a starting point for a 
phenomenological theory of "entanglement dynamics." Another variable of interest might 
be the repton density, which measures length per unit length [^. It may be possible to 



develop a hydrodynamics for this field as well [18 
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Figure Captions 

Fig. 1. Velocity fields and nematic directors witti various gradients and molecular orien- 
tations, a) and c) Gradient in the 2;-direction oiv±. This will rotate the director 
for all values of S. b) and d) Gradient in the x-direction of Vz- Only in d) will 
this rotate the molecule, i.e. when 5" < 1. 



16 



a) 



b) 



A 



A 



/\ 



c) 




d) 



/\ 



A 



Figure 1 



